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1.  Syanaij 


This  pap«r  deals  with  time  series  which  are  stationary  in  the  mean  but 
which  have  a  nonstatlonary  convarlance  function  which  is  described  by  a 
certain  parametric  SKxlel.  The  non-stationarity  results  from  ohanges  over  time 
of  the  parameters  In  the  isodel. 

The  first  part  of  this  paper  diseusses  the  concepts  of  stationarltyf 
covarisnce,  and  the  speotrm.  The  concept  of  temporal  spectrum  is  introduced 
and  applied  to  the  analysis  of  the  nonetationaxy  time  series. 

The  second  part  deals  with  the  parametric  model*  Statistics  related 
to  the  rate  of  sero~axls  crossings  and  the  rate  of  relative  maxima  and  minima 
of  the  time  series  are  used  to  estimate  the  parameters  in  the  model*  These 
estisiates  are  easy  to  compute  from  discrete  observations  of  the  time  series. 
Expressions  for  the  variances  of  these  eetljaators  are  given  ae  a  fxmotion  of 
the  parameters  and  the  sampling  period. 

2,  The  Spectrum 

In  the  analysis  of  time  series  data,  a  major  problem  is  the  character- 
isatlon  of  the  time  series  by  means  of  some  parameters  or  functions.  For 
example,  if  X(t)  denotes  the  signal  picked  up  by  the  antenna  of  an  FM  radio, 
the  AGC  (automatio  gain  control )circuit  estimates  the  average  power  E(X^(t)) 
and  corrects  the  gain  to  equalize  variations  over  time  in  this  parameter.  In 
order  for  this  system  to  be  effective,  the  average  power  parameter  must  be 
varying  slowly  enough  so  that  we  can  obtain  reasonably  accurate  estimates 
of  the  parameter  given  a  finite  record  of  X(t).  We  must  be  able  to  distinguish 
ohanges  in  the  characterization  parameters  from  transient  changes,  just 
"ordinary  noise, "of  the  time  series. 
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The  quest  for  slowly  varying  basic  parametere  leads  us  to  the  concept 
of  stationarity.  A  time  series  X(t)  is  called  strictly  stationary  if  for  all 
n,  the  Joint  distribution  of  {X(t^),  ...,X(t^))  is  the  sane  as  that  for 
{X(t^  ♦  h),,,.,X(t^  ♦  h)}  for  all  h,  Doob  [2],  X(t)  is  called  weakly 
stationary  if  the  density  of  X(t)  is  the  same  for  all  t,  and  for  all  t^  and 
the  Joint  density  of  (X(t^)f  X(t2))depends  only  on  |  t^  -  t^  I  •  Then  the  power 
spectrum  S(f)  and  the  autocovaziance  function  R(7r )  are  functions  associated 
with  X(t),  where 

S(f)  - 

(1) 


’>30 


exp(2^7rftjR('t')dr 


R(r)  -  E  (x(t+r)  x(t)) 


With  the  asauir.ption  that 

R(r) 


Oo 


d  t'<  00,  we  then  have 


2  /  cos  2:^f  2rs(f)df 
Jo 

and  by  differentiating  and  setting  have 


R"(0)  -  -(27/")^ 


(2) 


R^^^(o)  -  {zrrt  /f^s(f)df 

If  we  normalize  X(t)  such  that  E(X(t))  •  0  and  E(X^(t))  ■  1,  then  the 
even  function,  R(T')  ■  R(-7r)»  is  called  the  autocorrelation  function  and 
S(f)  is  called  the  power  spectral  density.  The  word  density  applies  since 
/  S(f)df  -  E(X^(t))  from  (l)  and  (2.  and  thus  j  S(f)df  -  1, 

Given  a  finite  record  of  X(t),  of  length  T  seconds,  there  are  several  ways 
to  obtain  an  estimate  of  the  spectrum,  see  Bartlett  [1]  and  Grenander  and 
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Rosenblatt  [3]«  Thsso  digital  prooodurss  or  a  good  analog  speotral  analgrasr 
glYss,  as  a  function  of  the  observations  In  the  record  of  X(t),  a  randoa 
function  S,j,(f)  such  that  for  aiqr  band  (f^^  f^),  (Middleton  [I4]) 


(3) 


S^(f)df  S(f)df  as  T  —  00  . 

f/v 


It  can  be  ahovn  that  If  we  pass  X(t)  through  a  linear,  tine  Invariant 

filter  which  filters  out  (analogous  to  an  optical  filter)  all  frequencies 

except  those  In  the  band  (f^f.),  then  with  X.  (t)  denoting  this 

tr  1  ^0  1 


filtered  tine  aeries, 


(4) 


S(f)df  -  E(I?  ,  :t)  )  . 


The  expected  \’alus  of  the  square  cf  a  tine  series  Is  called  the  average  power 
of  that  series.  Thus  (4)  shows  that  the  area  under  S(f)  for  f  In  (f^,  f^} 
gives  the  average  power  of  X(t)  In  the  band,  and  (3)  shows  that  the  Integral 
of  S^(f)  over  Is  a  consistent  estlaats  (WllJcs  [13]  )  of  this  average 

power* 

Moreover,  If  X(t)  Is  ergodlo  (Doob  [  2]  }  then  with  probabllit/  one 


(5) 


®{X(t)}  y*  X(t)dt 

T-r 

R  (r)  -  r  X(t  ♦r)  X  (t)  dt 
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B]r  using  the  Chsbjohsv  insquality  and  ths  Csussisn  nsturs  of  X(t)f  w«  hsvo 
froi  (5) 

T-r 


(6)  iL(r)  -T'V  x(t*r)x(t)dt 
-R(r)^Op  {rh 


whsrs  R,j,  is  called  the  sample  autocorrelation  of  X(t)  and  Op(f"^)  0 

in  probability.  Thus  from  (4)  and  (6)  we  have 


(7)  [At  •  f  *  On 

Jo  "o-'l  \  " 

which  suggeats  that  we  use  the  above  tins  average  of  ^  ,  ths  filtered 


time  series,  as  an  estimator  S^(f}  of  the  average  power  for  f  in  the  band 
(^0*  ^1^* 

To  do  spectral  estimation  as  above,  using  analog  msthods,  requires  a 
set  of  narrowband  filters  covering  a  lai^e  frequency  range. 

But  suppose  that  X(t)  is  stationary  in  the  mean  but  not  in  the  covariance, 
or  more  specifically  that  although  E(X(t))  •  0  for  all  t,  the  joint  density 
of  (X(t),  X(t  *^))  depends  on  t.  In  this  case,  S,j,(f)  is  an  estimate  of  ths 
average  power  for  f  is  a  band,  but  aimraged  over  distribution  changes  during 
the  sampling  period.  The  basic  diffioulty  is  that  we  must  distinguish  variations 
in  X(t)  due  to^ordlnary  noise”  versus  changes  over  tims  in  the  stnioture  of  ths 
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timm  ■•rlM.  In  order  to  handle  the  prohLea  we  would  like  to  ohareoteriee 
the  nonetatlonar/  X(t)  tagr  eone  generaliaatlon  of  the  power  epeetruBf  aueh 
ae  propoeed  by  Silreman  [?]•  Other  exteneiona  of  epeetral  aethode  to 
nonetationarj  proceeeee  are  given  by  Page  [11]  and  Laapard  [12],  However^ 
we  will  treat  the  problem  in  eomeidiat  of  a  heuriatio  manner. 

We  will  reatriet  ouraelvea  to  atationarjr  Oauaaian  time  aeriea  X(t) 
with  S(X(t))  ■  0  and  E(X^(t))  *  1  for  all  t^  and  auoh  that  given  a  record 
of  length  T  unita  the  joint  denaitj  of  (X(t^)»  XCt^))  ia  a  funotion  of 
I  t2  -  t^  I  for  t^  and  t2  in  the  T  unit  period,  i.e.  X(t)  ia  weaklj 
atationarx  duHng  that  period. 

furtheraore,  aaauae  that  there  exiata  a  r*<ifor 
there  ia  a  conatant  A  and  an  integer  a  auoh  that 


(8) 


E(X(t2)X(t^))  <A 


^2- 


Thia  aaana  that  for  laga  greater  than  'TT^he  tine  aeriea  ia  approxinatelj 
unoorrelated . 

Now  define  a  function  S(f),  called  the  temporal  apectral  denait/,  auoh 
that  aa  a  generalisation  of  the  definitlona  (1)  and  (2) 


(9) 


S(f) 


2/  coa  2  7rfrR('ir)d 

'V 


where  for  |'t^|  <  T 

(10)  R(r)  -  E(X(t  +  'V  )x(t)) 
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It  can  bt  shown  using  (7)«  (6)«  (9)«&nd  (10)  that 


(11) 


2 

S(f)df  -  E(X.  . 


(t))  ♦  0(T‘^) 


•/  t 

Jg  0*^  » 


whsrs  X.  .  Is  ths  rssult  of  s  linear  narrowband  filtering  operation  on 

V^l 

X(t)  for  the  (f^,  f^^)  band. 


MoreoTsr  for 

(12)  R(^') 


I't'l 

A  2 


* 

<  'X  bjr  taking  the  invsree  eoeine 


2tr  f  rs(f)df 


tranefonif 


for  the  tsnporal  spectrua  S(f),  and  given  the  condition  that  E(X^(t))  • 
we  then  have 


-00 


.00 

S(f)df  A  1 


e 


Suppose  we  have  a  finite  reoord  of  X(t)  during  which  X(t)  is  etationarj 
and  we  wish  to  estinate  the  temporal  epeotral  densitgr  S(f)  for  f  in  sobs 
frequencj  band,  ’’e  could  use  the  method  of  narrowband  filtering  (see  (?)  and 
(11)  )  to  estimate  S(f).  However,  in  nanj  applications  the  period  of  stationarj 
is  so  short  as  to  make  the  spectral  estimates  which  are  fbnotions  of  observe* 
tions  taken  in  the  period,  quite  inaccurate.  For  example,  the  ambient  sea 
noise,  as  measured  by  hydrophones  placed  in  the  deep  ocean,  is  often  quite 
nonstationary  due  to  sudden  changes  in  the  paths  of  propagation  of 
sounds  in  the  ocean.  If  %fa  wish 
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Acoiirate  spectral  estimates  for  f  In  a  relatively  wide  band«  in  many  applica¬ 
tions  we  have  to  take  a  record  of  X(t)  which  is  longer  than  the  period  of 
stationarity,  and  thus  the  temporal  spectnn  approach  is  not  meaningful  since 
we  have  averaged  over  a  change  in  structure. 

However^  suppose  we  know  from  theoretical  consideratims,  or  fras  past 
sxperimentation,  that  when  the  time  series  has  a  period  of  statlonarity, 
the  temporal  spectral  density  S(f)  is  one  of  a  certain  parametric  family  of 
density  functions.  For  example,  suppose  we  can  say  that 

(13)  3(f).  J^T) 

where  r  and  X  are  two  unknam  parameters.  In  this  case  we  have  S(f) 
described  by  a  family  of  Gamma  densities  (Figures  1  and  2).  For  different 
stationarlty  perixls,  r  and  X  take  different  values.  We  thus  have  reduced 
the  nonstationarlty  of  X(t)  to  the  changes  of  two  (or  in  general  a  finite 
set)  parameters  in  the  temporal  spectrum.  From  (12)  we  see  that  the  auto¬ 
correlation  function  R('tr)  la  a  nonstationary  function  of  these  time  varying 
parameters. 

However,  a  stationary  X(t)  with  a  spectral  density  as  in  (13)  is 

deterministic  (Doob  [2])  since 

00 

(1  ♦  f^)”^  log  S(f)df  -  -  00 

*  V 

This  means  that  riven  a  time  t  ,  the  random  variable  X(t  )  can  be  estimated 

* 

with  zero  variance  by  a  random  variable  Y(t  )  which  is  the  result  of  a  certain 
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linear  operation  on  the  record  of  the  infinite  paetf  i.e,  on  X(t)  for  all 

* 

t  <  t  .  But  el  nee  we  are  dealing  with  nonetationarjr  time  aeries^  reaults 
baeed  upon  obeenratlons  over  the  infinite  past  do  not  eeeei  to  be  relevant* 

3*  Parametric  Model  for  the  Temporal  Spectxnna 

Suppose  we  make  dlacrete  observations  of  a  Oauaelan  time  eeriea  X(t) 
dxiring  a  stationary  period  where  the  temporal  apeotrum  S(f}  is  a  Gamma  density 
parameterized  by  r  and  \  -  ae  is  given  by  (13)#  l.e. 

(13)*  S(f)  -  ^  ° 

We  wish  to  estimate  r  and  k  from  the  discrete  ''bservatione*  The  eatlmators 

are  of  no  use  if^  in  order  to  obtain  relatively  accurate  estimates  of  r  and 

we  need  a  sampling  period  longer  than  the  period  of  stationarity  of  X*  Any 

estimators  r  and  X  based  on  observations  taken  over  time  from  X(t)  will 

have  two  components  of  randomness,  the  first  due  to  the  stationary  randem 

fluctuations  over  time  of  X(t),  and  the  second  due  to  the  changes  in  the 

covariance  of  X(t)  vdiich  recur  from  time-to-time. 

The  rate  of  zero-axis  crossings  of  X(t)  in  an  interval  (0,T)  is  a 

consistent  estimator  of  3  ■  2[^  f^S(f)df]  as  T-»od,  and  the  rate 

of  relative  maxiina  and  minima  of  X  in  (0,T)  is  a  consistent  estimator  of 

•  2[  J  f^S(f)df/  J  f^S(f)df  .  We  will  express  r  and  X  as  functions 

*  *  * 

of  3  and  3  ,  r  ■  f(3,3  )  and  X  ■  g(3,3  )»  and  estimate  r  and  X  by 
r  •  f(3»3  )  and  X  -  g(3,3  )  where  3  le  related  to  the  rate  of  aero-crosoir^s 

/V  * 

and  3  is  related  to  the  rate  of  maxima  and  minima,  but  they  are  easy  to 
compute  from  discrete  clipping  of  X(t),  which  we  will  describe  later*  The 
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consistency  of  the  estimators  Is  not  important  due  to  the  nonstatlonarity  sf 

X(t),  For  the  estimators  to  be  of  real  value  in  this  context,  their  variances 

should  decrease  relatively  rapidly  with  increase  in  the  period  of  sampling. 

* 

In  this  Mork  we  will  obtain  the  variances  for  the  estimators  of  3  and  3  as 

functions  of  r,  X,  and  T.  For  the  case  of  low  frequencies,  that  is  for  the 

mean  frequency  J*  | f  f  S(f)df  -  r/X  in  the  range  20  cpe  to  UO  cps, 

calculations  (Hinich  C  /»])  show  that  the  standard  deviation  of  the  estimates 
* 

of  3  And  3  is  less  than  10^  of  3  for  a  sampling  period  of  two  seconds  and 

less  than  of  3  for  a  sampling  period  of  ten  seconds. 

Suppose  we  sample  X(t)  in  a  discrete  manner  by  observing  every  tT  time 

units  whether  the  time  series  is  positive  or  negative.  If  X  is  positive, 

we  mark  a  1,  if  negative  we  mark  a  -  1,  i.e.  we  define  a  discrete  time 
« 

series  X  where  for  each  integer  k 

,1  if  XCkT)  >  0 

(U)  X*(kr)-  J 

^-1  if  X(k^)  <  0 

This  is  simply  the  discrete  version  of  infinite  clipping  (Lawson  and  Uhlen- 
beck  [  5]  ).  ^ith  the  advent  of  high  speed  sampling  techniques  using  digital 

electronics,  it  is  a  standard  procedure  to  subject  a  random  process  to 

infinite  clipping  in  order  to  reduce  the  data  and  put  it  into  a  convenient 
form  (a  string  of  binary  numbers)  for  real-time  analysis.  The  estimators 
for  this  parametric  model  have  the  useful  property  that  they  are  easy  to 
obtain  fran  the  clipped  process. 
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Now  let  li{T  denote  the  number  of  tranaltione  of  X  during  the 
Interval  (0,T)  from  -1  to  1  and  from  1  to  -1.  Fron  Leona  3  we  have 

(15)  00  Nd.-V  )A  -  P  ♦  0  i't^) 

00 

where  -  4  |  f^  S(f)df  and  thus  we  have  fron  (13), 

J  -00 

2  —2 

3  "LX  r(r*l).  Thus  from  (15),  we  have 

(16)  li«T-*oo  N(T,r)A  - Cr(r  ♦  1)]  ♦  0  (-tr^) 

* 

The  rate  of  binary  traneitions  of  X  -  N(TjT')A  “  i®  intimately 
conncected  with  the  rate  of  eero-axia  oroeeings  of  X  provided  the  proceaa 
X  is  'well-behaved.'  To  argue  heuriatically,  select  'X'  sufficiently  aaall 
such  that  in  any  'tf  unit  interval,  the  probability  of  two  or  more  axis- 
crossings  can  be  neglected.  Since  either  there  is  one  crossing  or  no 
crossing  in  each  interval  {k  't' ,k  't'  number  of  axis-crossings 

of  X  in  the  broad  interval  (0,T)  is  Just  N(T^'V)  and  thus  N(T^'l:'  )A  i®  Lhe 
rate  of  axis-crossings  of  the  time  series* 

In  ordor  to  give  the  variance  of  N(T,'^)A  we  need  the  autocorrelation 
R(T)  of  X(t)  .  By  taking  the  cosine  transformation  of  (13),  we  have 

(17)  R(Tr)  -  Re  (1  - 

where  i  is  the  complex  unit  and  Re  z  is  the  real  part  of  a  complex  a.  We 
shall  restrict  r  to  be  an  integer  for  computational  convenience  in  obtaining 
the  sampling  variances  of  several  estimators. 


Using  (17)  in  Lemna  U  along  with  the  relations 


(18) 

where  I 
second 

(19) 


where 

(20«) 

(2a>) 


r’  (TT)  --2Trl  -('*!) 

R*  (-f-)  .  .  (2  tT)^  -  Re  (1  -  ^-^3:  )  -  (r»2) 

I  ^ 

i  is  the  first  derirative  of  R  with  respect  to  o  snd  R"  is  the 
derivative^  and  changing  variables  t  •  ZW^x,  we  have 


Var  [N(T  , -r  )A] 

StTT* 


■  2^  f  I- 


(1  -  2^  )  .,(t)dt  *  (^) 


T*  .  X  X 


B-f  [r(r*l))^/* 


•Jt)  e 


Cb  (t)] 


1/2 


1  -  (t) 


j  [1  ♦  c  (t)arctan  c  (t))  -  1 

^  t^\  r  r 


2  2 

bj.(t)  -  [1  -  e^  it)  ]  [1  -  e^2 
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(20c) 


_r__  2  2 

r+1  (t)  - 

[l-«!  (t)l^^^  Cb,(t))^/^ 

r  r 


(20d) 


d.(t)  ■  r“^  (rx)  - 


r*  r*l  ,  .  2J41 

I  (,.*,)  (-i)^(^) 

1-0 


[1  ♦  (7)*  j 


(20e) 


•^(t)  -  R(rx) 
r 


j.* 

£  (y.)  (-ly  (;) 

a-jlJL—j: _ 


tl  ♦  (;)  *)  ' 


wher« 

* 

r  • 


r 

2 


if  r  is  ev«n 


2 


if 


r  is  odd 


Renembsr  that  rA 


S(f)df  which  is  the  "m»an"  frsqusno/  of  ths 


tenporal  spectrum  , 

We  will  now  discuss  the  estimator  related  to  the  number  of  maxima  and 
minima  per  unit  interval.  Define  the  derivative  process 
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(a) 


X*(t) 


X(t»h)  -  I(t) 
h 


where  the  oonwergence  is  In  the  mean  (Doob  C  2]  )•  If  X  le  a  Gaussian 
process  with  an  autocorrelation  R  such  that 

(22)  RC-r)  -  1  ♦  R^CO)  ^  ♦  R^^^ 

almost  everywhere  and  Is  also  Gaussian  (Doob  [  2]  ).  Moreover,  it  can  be 

I  *»  •» 

easily  shown  that  the  autocovarlar.':e  of  X  is  -R  where  R  is  the 

I 

second  derivative  of  R  with  respect  to  L  •  If  we  normalise  X  by 
defining  a  process 

(23)  y(t)  -  [-r"(o)  x'  (t) 

ws  have  E{T(t)}  •  0,  E{Y^(t)}  ■  1  and 

(24)  E{Y(t4  X  )Y(t)}  - 

R  (0) 

If  we  let  D(f)  denote  the  spectral  density  of  Y,  we  have  from  (2)  and  (24) 


(25)  ^A(f)df  -  (2Yr)"2  [R'*(o)r^  R^^^  (0) 

-  /f^  S(f)df/  /f^S(f)df 
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The  process  X  has  a  relative  maxima  or  minima  at  t  if  and  onlj  if 
X  (t)  ■  0,  and  thus  If  and  only  if  y(t)  •  0,  If  we  let  M(T^  be  the 
number  of  transitions  from  1  to  >1  and  from  -1  to  1  of  the  diaorete 
process 

^1  if  Y(k'r)  >  0 

(26)  Y^lkT)  - 

(-1  if  Y  (k'T)  <  0 

for  k  ■  0,1, ,  then  for  small  T*  ,  M  is  approximately  the  number 
of  relative  maxima  or  minima  in  the  interval  (0,T).  Therefore  M(T,  X 
is  essentially  the  rate  of  maxima  and  minima  of  X.  Using  (23)#  (25)»  and 

(26)  in  Lenina  3*  we  have 

(27)  -  3*  ♦  0  (r^) 

-  2[Xf^s(f)df/Xf^s(f)df]^/^  ♦  o(r^) 

* 

where  3  is  ^’be  expected  niinber  of  relative  maxima  or  minima  per  unit 
Interval  (Rice  [  7]  ). 

a 

However  we  do  not  have  to  differentiate  X  in  order  to  estimate  3  •  Let 
us  define  another  discrete  process 

^  1  if  X((k^l)  Tr  )  >  X(kr) 

(28)  X**(kr).  j 

(^-1  ifx((k>i)‘r)  <x(kr) 

Fran  the  definition  of  X  and  (28),  it  is  fairly  clear  that  the  rate  of 
binary  transitions  of  X  for  sufficiently  small  is  also  M(T  j't')/T. 
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Applying  (27)  to  the  c&se  of  e  Ganme  speotrunif  we  heve 

2 

(29)  [(r*2)  (r^3)  ) 
since  from  (13) 

yW(f)<lf  -  X''‘r(r*l)(r*2)(r*3) 

We  can  obtain  the  Yarlance  of  M(T,'C')/r  fairly  easily  based  on  the  deriYS- 
tion  of  (19).  Applying  (16)  to  {2U),  we  have  for  the  autooorrelation  of  T, 

(30)  E{Y(t^'r)y(t))  -  Red-  — ^ 

which,  from  (17),  is  Jxist  the  autocorrelation  of  X  with  r  replaced  by  r*2. 
Setting  —  T**  -  T  and  3*  -  2\’’^  [(r^2)(r*3)}^^^  (just  3  with 

r*2  instead  of  r),  we  hare 


(31) 


lim  Var[M(T,  T  )/T] 
X-¥0 


(t)dt 


‘  (^)  3 

which  is  the  variance  of  N(T,  'Z'  )/T  as  given  by  (19)  and  (20),  but  whore 
r  is  repr'*''ed  by  r+2. 
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We  will  now  propose  estimAtors  or  r  and  From  (16)  and  (29)  we  have 


which  jlelds 


2  0^  *2  « 
(32)  r-^ _ iJiS.  U4a_El 

2 

2(0  -  n 

Let  the  estimator  of  r  by  r  where 


A  ,  SN^-M^  ♦  [M^  ♦  IL 
2(M^-N^) 

for  N  •  N(T^'lO  and  M  -  M(T  ,T).  Thus  from  (16)  and  (29)  we  see  that  r  ie 
a  ooneletent  estimator  of  r.  Since  r  le  reatrioted  to  be  an  integer^  ohooee 
[r]  as  the  estimate  of  r  where  [r]  Is  the  closest  integer  to  r.  [r]  la  alao 
a  consistent  estimator  of  r  given  the  statlonarlty  of  X,  But  einoe  X  is  in 
reality  nonstationary,  the  consistency  of  the  estimator  Is  indeed  a  weak 
property.  The  usefulness  of  r  will  be  mainly  determined  by  the  Tarianoe  of 
r,  which  we  do  not  give  since  we  have  not  developed  the  covariance  of  N  and 
M.  However,  if  in  a  certain  application  for  a  specified  pair  of  parameters 
r  and  X,  the  variances  of  N/T  and  M/T  are  reasonably  small,  then  it  is 


worthwhile  to  find  the  exact  variance  of  r  as  a  function  of  T  by  Monte  Carlo 
techniques.  In  the  case  of  mean  frequencies  20  cpe  <  rA  ope,  [r]  had  a  Sf* 
standard  deviation  for  T  of  less  than  a  minute. 


Again  from  (16), 


"-fer 

la  a  cona latent  eatlmator  of 

/N 

Both  r  and  \  are  relatively  easy  to  compute  from  a  eequanoa  of  dlaoreta 
obaervatlona  on  the  clipped  version  of  the  time  aeries  X(t),  Slnca  in  a  grsat 
■any  applications  dealing  with  the  above  type  of  nonatatlonary  time  aeries, 
the  process  Is  first  subjected  to  Infinite  clipping,  the  estinatlon  BMthod 
given  above  offers  a  practical  procedure  tar  dealing  with  the  changes  of  tha 
hamonj.e  power  levels  over  time. 

4«  Eraodic  Reeulte  for  the  Rate  of  Axle-Croeainga 

Let  N(T)  be  the  number  of  tero^axls  crossings  of  X(t)  in  the  interval  (0,T) 
and  aaaumlng  that  E{X(t)}  ■  0  and  E{n(T))  <  00,  define 

(34)  I.  . 

At  first  glancs  It  looks  as  if  3  is  a  function  of  T,  If  X(t)  Is  strictly 
stationary  than  for  all  T,  the  expected  number  of  sertHcrossings  in  any 
T  unit  interval  Is  the  same  as  £{N(T)).  Now  if  we  choose  any  unit  of  time 
‘t’  and  let  K  be  the  positive  integer  K  ■  T/t',  then  from  (3A) 
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and  thus  0  is  not  a  function  of  T. 

Rica  [7]  shows  that  if  X  is  Gaussian^  0  •  2[  >^f^3(f )df td^art  S(f) 
is  the  spectnuB  of  X.  However,  in  Rioe*s  derivation  and  also  In  MoFadden's  [6], 
it  is  assuBied  that  a  X  can  be  found  suffioientlj  small  such  that  in  the 
interval  (t,  t  *  't),  the  probability  of  exactly  one  axis-crossing  is  simply  0 
and  the  probability  of  two  or  more  crossings  can  be  nsglsctsd*  This  assumption 
should  be  carefully  considered  sines  if  X  is  white  noise,  there  ars  infinitely 
many  seres  in  any  interval.  We  will  now  give  three  simple  Isiss  which  will 
show  that  N(T)/T  is  a  consistent  estimate  of  0;  that  the  second  moatent  of  the 
speetnim  of  a  Gaussian  process  (if  it  exists)  can  be  estimated  from  diserete 
clippingj  and  then  give  the  variance  of  the  consistent  estimator. 

Lenaa  li  Let  X(t)  be  a  strictly  stationary  and  metrically  transitive  (Doob[2]) 
random  process  and  N(T)  be  the  number  of  axis-crossings  in  the  interval  (0,T). 
Assume  that  E{n(T))  exists.  Then  as  T  — »  co 

N(T)A-8 

with  probability  one. 

Proof  I  Express  T  in  integer  multiples  of  a  convenient  unit'T,i.e.  T  ■  for 

a  positive  integer  K.  Let  denote  the  number  of  axis-crossings  of  X  in  the 

interval  ((k-l)'l'  ,k  'V)  for  k«l,,..,K  Thus  we  have 

K 

(35)  N(T)  -  \ 

Since  X  is  strictly  stationary  and  metrically  transitive,  so  is  ths  diserete 
process  Thus  we  have  the  ergodio  result  that  as  K-»oo,  with  probability 

ene 
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(36) 


k-1 


But  fro«  (34),  “  P  "C"  4nd  thus  by  using  (35)  in  (36)  and  diYiding 

by  ^  we  have  the  desired  result. 

Lenana  2t  Let  X(t)  be  a  strictly  stationary  and  metrically  transitive 
random  process.  Given  a  T  ,  define  the  discrete  process 

•1  if  X(k  T)  >  0 


X  (k  'T  )  - 


-1  if  X(k  r )  <  0 


for  the  integers  k  ■  0.  1,  ...  .  Let  N(Ty  'IT  )  denote  the  number  of  transi- 
tions  of  X  from  1  to  -1  and  from  -1  to  1  in  the  interval  (0,T)  where 
T  •  K  t'  .  Let  R  (  'L")  •  E{X  (k 't' )  X  ((k^l)  TT)}  which  does  not  depend 
on  k  since  X  is  strictly  stationary  and  metrically  transitive*  Then 


N(T,  t)/t-  (2  r)  [1  -  R*('r )] 


with  probability  one  as  T  os . 

* 

Proof:  By  the  ergodicity  of  X  ,  the  sample  autocorrelation 


K— 1 

(37)  R*(ir)~^  L  X*  ((k-1)  T)  X*  (k  r) 

T  k-0 


converges  with  probability  one  to  the  autocorrelation  R  asK«T/'^— oo* 

0 

Applying  the  definition  of  X  to  (37)  we  have 
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(  XT)  -  K‘^  [K-2N(T,  'C)  ] 

and  thus  bj  dividing  bj  ^  in  the  convergence  R.^,  -»  R  ,  we  have  the  deaired 
result. 

Now  suppose  X(t)  is  stationary  Gaussian  process  where  E{X(t))  ■  0  and 

2 

(t))  ■  1.  It  then  possesses  the  properties  required  for  the  If  ss. 
Thus  N(T)/r  3  with  probability  one  as  T-»  oe  and  siailarlj'  for  laiBa  2. 
Slepian  [6]  we  can  show  Rice's  [7]  result  that 

(3,)  e-Tf  [-r-(0)]V2  -  2  C  i 

given  that  R  is  t%d.ce  continuously  differentiable  at  ^  ■  o. 

But  lawson  and  Uhlenbeck  C$]  show  that 

(W)  h'(  -T)  -  2  MMin  R(  r  ) 

Bjr  .erl*.  ucpuisicn  of  .rooln  In  poworo  of  ,  wo  out  ohow  fro.  (40) 

(U)  ^  (-R"(0))^''* 

»  (R-CO)]^  -  R^'*^  (0) 

[-R*(0))^'* 

where  now  we  assuae  that  R  has  the  expansion 
(42)  R(tr  )  -  1  ♦  R"(0)  (0)  ^ 


24  "ff 


0  {'th 


Trm 


from  (39)«  (41)  and  Lesna  2  we  hav<^ 
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La—  3t  Suppose  X(t)  is  s  stationary  Oaussian  procsss  with  an  auto- 
oorralatien  function  R(  T' )  as  given  by  (42) •  Than  with  probability  ona 

^-0  '*(T.r)A-P 


and 

,  2 

[HCr.rjA-P)  .yiScDdf} 

»  0  (-r^) 

whara  S(f)  is  the  spectral  density  of  X. 

Wa  will  not  derive  the  variance  of  N(T,  )/7  as  a  function  of  T, 
and  3*  From  (33),  Lemma  3,  and  K  •  T/'^  we  have 


(43)  ECRj(Tr)]^-l 


43  r  ♦  4  E[N(T, 'jr  )A]^ 
♦  0  ('T^) 


anl  from  (37) 


(44) 


E[R^  ('T)]^ 


K—1 

E  E{x*(k^i)  r)  X*(k  Tr)x*((j^i)  r  )X*(J  -jr)) 
J-i 


Let  D(x,  't')  denote  the  conditional  probability  that  X(t)  erosaas  the 
axis  in  the  interval  (x,x^  't' )  given  that  it  crossed  in  (0,^).  For 
Gaussian  X(t),  Rice  [?]  derlvee  a  function  U(x)  such  that 

c 


D(x,  ^)  •  u(x)  'ir  *  0{  'Ir^)  and 


U(x)dx  -  1  as  e  -►  0, 
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Sine*  the  probability  of  the  procese  croselng  the  axle  In  the  intenral 
(0,  Tr)  we  then  have  for  |  k  -  J  j  ■  n  >  1 

U5a)  E{x’‘(k*l)  'r)X*(k  ^)X*((J^1)  'tr)X*(j  ir)) 

-  U(n  'T)  't')  -  2[1-U(n  ^)Tr(P'2^)] 

♦  {l-U(n-r)  (BTr)  -2[l-u(n 

♦  Qi^cr^) 

-  1  -  43  Tr-  ♦  4  u(n  Ir)^  ♦  0(  TT^) 
and  for  k  ■  J 

(45b)  E{X*(k>l)  r')X*(k  'Zr)X*(  )X*(j  Tt)  -  1 

Setting  (45)  in  (44)  and  eunaing,  we  have 


(46)  E(R^(T')]^  -  1-437:  ♦  60  T"^  ^i^(l  -  ^)  U(n7r)7r 


n«l 


♦  0( 


and  thua  fron  (43) 


(47)  ECN(T,  t  )A]*  - -I®-  E  (1 )u(n't)'t' 

n-1 

Letting  '7"  go  to  aero  In  (47)  and  since  11«^  ^ 


U(x)dx 


1 

2 


we  have 

(48)  Var[N(T,  r)A) 


'’T 

_2£L 

T 

w/  0» 


(1-  f )  U(x)dx 


♦ 
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If  X  !•  OA\i0ti«n  with  AutoeorrelAtion  R(7r')*  using  Rlcs's  derivation  of  U(x), 
we  have 

Lcnnna  Ui 

0  Var[N(T,  Tr  )A]  -  I  (l-f)A(x)[l  ♦  D(x)arctan  a(x)>lx 

Jo 

71  , 

where 

*(x)  -  tB*(x)  -  C*(x)]  [1-  R*(x)r^/* 

B(x)  -  1  -  R*(x)  ♦  [R"(0)r^  [R'(x)J* 

C(x)  •  R(x)  tR’(0)r^  tR'(x)]^  ♦  R"(x)  [R’CO)]*^  tl-R*(x)) 

D(x)  -  C(x)  [B*(x)  -  C*(x)]  • 


This  Is  slallar  to  the  result  of  Steinberg,  et.al.  [10] 


300 


1 - 1  I  I - 
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FIG.  2  GAMMA  SPECTRUM  FOR  r«5 
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